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We derive a semiclassical scheme for the conductance through a rectangular cavity. The trans- 
mission amplitudes are expressed as a sum over families of trajectories rather than a sum over 
isolated trajectories. The contributing families are obtained from the evaluation of a finite number 
of continued fractions. We find that, contrary to the chaotic case, the conductance fluctuations 
increase with the incoming energy and the correlation function exhibits a singularity at the origin. 
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I. INTRODUCTION 

The density of states of small quantum systems is 
known to be sensitive to the underlying classical mechan- 
ics. The chaotic or integrable character of the classi- 
cal dynamics translates into the statistical properties of 
the spectra 0. The spectral correlations of classically 
chaotic systems are described according to random ma- 
trix theory Q and the fluctuations with respect to the 
mean density of states are less pronounced than for in- 
tegrable systems. The connection between classical and 
quantum properties is particularly clear within the semi- 
classical approximation: the density of states of classi- 
cally chaotic systems is expressed in the Gutzwiller trace 
formula || as a sum over [isolated] periodic orbits, while 
integrable systems are described by the Berry- Tabor for- 
mula Q , where the the density of states is given as a sum 
over families of trajectories (or invariant tori). 

What is the corresponding situation in open systems? 
Transport properties instead of spectral ones are usually 
considered. The transmission amplitudes and probabil- 
ities are given as semiclassical expansions over trajecto- 
ries that traverse the scattering region providing 
a link between classical and quantum properties. Such 
an approach is firmly established in the case of isolated 
trajectories, as typically found for chaotic dynamics, but 
only recently the subtleties associated with an integrable 
geometry have started to be addressed J?] |llj . As in 
the Berry- Tabor formula for the density of states, the 
existence of families of trajectories degenerate in action 
modifies the structure of the semiclassical expansion. In- 
tegrable cavities fed by leads present the additional prob- 
lem that conserved quantities in the scattering region are 
not necessarily conserved in the leads. The semiclassi- 
cal description of transport through integrable cavities is 
a necessary step towards understanding the influence of 
the underlying classical mechanics in open systems. This 
is the goal of this work, where we present a semiclassi- 
cal evaluation of the transmission amplitude through a 
rectangular cavity and apply it to the study of the con- 
ductance fluctuations. 

The interest in open systems stems from their parallel 



with closed systems and also from the physically mea- 
surable realizations that have been developed in the last 
few years. High mobility-semiconductor microstructures 
at low temperatures give access to the ballistic regime, 
where transport is dominated by the geometrical scat- 
tering of electrons off the walls of the lithographically 
defined cavities [jl2p|l4]]. On the other hand, microwave 
cavities provide an easily tunable system to study the 
propagation of electromagnetic waves in a defect-free re- 

gion g. 

Quantum interference in ballistic cavities gives rise to 
conductance fluctuations under a small perturbation of 
the system (magnetic field, Fermi energy of the incoming 
electrons, shape of the cavity, etc) |l^|l6|,[L7| . For chaotic 
cavities, the characteristic correlation lengths of the fluc- 
tuations are well described by the existing semiclassical 
theories while the universal character of the vari- 

ance has been addressed within a Dyson hypothesis for 
the scattering matrix |n],|2(]]. Random matrix theories 
are obviously not appropriate to treat integrable systems, 
forcing us to rely on semiclassical expansions and nu- 
merical calculations. The early experiments on ballistic 
cavities [[l2| yielded a more structured spectrum of con- 
ductance fluctuations for the integrable case than for the 
chaotic one. This tendency was reproduced by numerical 
calculations showing that in rectangular cavities the 
higher harmonics of the conductance spectrum are more 
pronounced than in the chaotic case. Also, the variance 
of the conductance appeared as increasing with the mode 
number. The lack of transport theories for integrable sys- 
tems (which do not enjoy the universality properties of 
chaotic cavities) and the difficulties with their fabrication 
are probably responsible for the fact that subsequent ef- 
forts on conductance fluctuations were mainly centered 
on the chaotic case [^l[ . 

Weak-localization, the decrease of the average conduc- 
tance due to constructive interference of time-reversed 
backscattering trajectories, is another quantum inter- 
ference effect which has recently been studied in bal- 
listic cavities. The underlying classical mechanics was 
proposed [g2| to translate into a different shape of 
the magnetoconductance peaks: a lorentzian line-shape 
for chaotic cavities and a triangular form for the in- 
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tegrable case. Chaotic cavities indeed exhibited the 
lorentzian line-shape while rectangular and circular ge- 
ometries showed triangular or lorentzian peaks depend- 
ing on the experiment ||l6| , |l7| , p3| - |2^] . The importance of 
the leads has been demonstrated pq] and the disagree- 
ments between different experiments are not settled yet. 
This illustrates the difficulties in dealing with integrable 
cavities and the importance of developing appropriate 
theoretical descriptions. 

In this work we consider a square cavity (of length 
side 1) with hard walls and connected to leads on oppo- 
site sides of the square by openings of size W (see inset 
Fig. . (For our purposes square cavities are no less gen- 
eral than rectangular ones since we can always scale one 
of the sides.) We develop in Sec. |l|a semiclassical expan- 
sion of the transmission amplitude between two modes in 
terms of families of trajectories degenerate in action. The 
determination of these families (Sec. Ill) results in a dio- 
phantic problem and the evaluation of a finite number 
of continued fractions. This constitutes a very efficient 
numerical method that is to be compared with the nu- 
merical quantum approaches based on discretization of 
real space or truncation of Hilbert space. The classical 
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conductance (Sec. tV) and the unbounded fluctuations 
(Sec. 0), increasing with the incoming energy, are ob- 
tained from the numerical implementation of the contin- 
ued fraction approach. The correlation function of the 
conductance, when properly normalized, is shown to be 
well defined in all the energy range and shows a cusp-like 
singularity at the origin. Further approximations on our 
semiclassical expansions allow us to demonstrate that the 
unbounded conductance fluctuations grow linearly with 
the incoming Fermi wave- vector k. In the conclusions 
(Sec. VI) we analyze the experimental relevance of our 
findings and the implications for other integrable geome- 
tries. 

Transport through a square cavity has recently been 
considered by Wirtz and collaborators ]Tl]| , where al- 
ternative forms of the transmission coefficient were pro- 
posed and the correspondence between quantum results 
and classical trajectories was established. The dynamics 
through rectangular cavities constitutes a paradigm for 
the interplay between scattering and integrability, and 
has been treated in other contexts. For instance, it has 
been presented by Zwanzig p7| as an example where the 
short memory approximation to the migration of a clas- 
sical dynamical system between regions of configuration 
space is entirely wrong. 



II. SEMICLASSICAL FORMULATION 

We consider a phase-coherent cavity connected to two 
reservoirs through leads that support N propagating 
modes. The conductance in the Landauer formalism is 
proportional to the transmission coefficient T 
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T : 
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(1) 
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The complex numbers tba are the transmission ampli- 
tudes from the mode a at the entrance lead to the mode 
b at the exit lead. They are related to the Green func- 
tion (evaluated at the Fermi energy E F ) connecting the 
point (0,y) at the entrance with (L,y') at the exiting 
lead through J28[: 



tba 



ih^V a V b 



dydy' 4>t(y') <Pa{y) 

*G(L,y';0,y;E F ) (2) 



where we have discarded an unimportant phase factor. 
4>a, b represent the transverse wave functions in the leads 
of width W {<j)a(y) = \f2/W sin(a7n//W)) and w Qjb the 
longitudinal velocities associated with the modes a and 
b. 

In the semiclassical approximation, the Green function 
is expressed as a sum over trajectories joining the points 
(0,y) and (L,y>) §: 



G scl (L,y / ;0,y;E F ) = 
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s(y,y') 

exp [jrS s (y', y, E F ) - i|/i s 



(3) 



S s is the classical action along the trajectory s. Since we 
are considering billiards without magnetic field S s — kL s , 
with k = mv/Ti the Fermi wavevector and L s the tra- 
jectory length. Denoting 9 and 6' the incoming and 
outgoing angles, the preexponential factor is given by 
D = (vcosd' /m)~ 1 \(d6/dy') y \. We include in the phase 
/i s the Maslov index counting the number of constant- 
energy conjugate points and the phase acquired at the 
bounces with the hard walls. 

The semiclassical expression of the transmission ampli- 
tude is obtained by using the above approximation of the 
Green function in Eq. (0) and performing the y and y' 
integrals to leading order in h. A stationary-phase inte- 
gration over y selects the trajectories entering the cavity 
with the angle 8 a (a — ±a and sin# a = an/kW) and 
leaving at the point y' , leading to 
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x exp 



iS s (y',6- a ,E F ) 



The reduced action is 

S(y', 6 a , E F ) = S{y\ y , E F ) + tnray /W 



(4) 



(5) 



where yo(0 a ,y') is the initial point selected from the 
stationary-phase integration. The prefactor is now given 
by B = (v cos9')~ 1 \(dy/dy')e\, and the Maslov index is 
increased by one if (d6/dy)y> is positive. 
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When nearby trajectories are non-degenerate in action, 
as is typically the case in chaotic cavities, the integral 
over y' can also be done by stationary-phase and the 
semiclassical t oa is expressed as a sum over trajectories 
with quantized initial and final angles j^J^] . On the other 
hand, when trajectories come in families with the same 
action (or length) the phase in Eq. (||) is linear in y', 
preventing another stationary-phase integration. This is 
the case of the trajectories directly crossing a ballistic 
cavity from the entrance to the exiting lead, where the 



y'-integral can be exactly done 



The trajectories 



going through a square cavity also come in families, or 
bundles JTlj ], and the integration over the j/'-coordinate 
is very similar to the case of direct trajectories since the 
dynamics in an extended space is that of free particles 
(Fig-0). 




FIG. 1. Unfolded space for the dynamics in a square cavity 
(inset). The trajectory entering at the lowest point (y = 0) of 
the left lead with an angle 6 is shown in the original and ex- 
tended space. It belongs to the family that leaves the cavity 
through the exiting lead (po, qo) and has a weight So. For the 
particular chosen 9 there are two families contributing to the 
transmission amplitude (0 and 1) and one family contributing 
to the reflection amplitude (0*). 

Symmetry arguments for a rectangular cavity dictate 
that tba — if a + b is odd. For even a + b we perform 
the y'-integration for each family n defined by the coor- 
dinates (p n , q n ) of the exiting lead in the extended space 
obtaining 
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x exp 
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y^ and y'^ are the extreme points of the exiting inter- 
val, e(q n ) =0 for even q n and s(q n ) = 1 for odd q n . The 



parity of q n appears in the phase due to mirror symme- 
tries involved in going to the extended space. The trajec- 
tories of the n-th family have a length L n = p n j cos 9 a (all 
lengths are expressed in units of the side of the square) . 
In the extended space we have free motion, therefore the 
phase jx n is simply given by the p n — 1 bounces with the 
vertical walls and the q n bounces with the horizontal 
ones; that is, fj, n = 2(p n + q n — l). The trajectories that 
contribute to the transmission amplitude are those going 
from the left to the right lead, therefore only the values n 
with odd p n should be considered. Conversely, the even 
values of p n yield the reflection amplitude. For transmis- 
sion amplitude we can simplify the phase of Eq. (^|) and 
write 



i I cos ( 

tba = — 777 



y^e n exp[ikL n ] 



(8) 
(9) 



W V cos(9 a 

x {I n (a+b) - I n (a-b)} 

We have defined the phase 

e n = exp[i7r(a + l)e((7„)] , 

and the reduced length 

L n = Pn cos 9 a + q n sin 6 a . (10) 

Similarly to the case of direct trajectories 0, a+b — 
has to be treated separately for the ^'-integration. How- 
ever, the corresponding result is included in Eq. (||) by 
taking the limit x — > 0. Obviously, a + b = corre- 
sponds to the maximum transmission since this is the 
case where the classical trajectory arrives to the exit- 
ing lead with the quantized angle of mode b. As em- 
phasized in Ref. the diffractive terms obtained for 
a + b =/= become more important for families with a 
weight 5„ 
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-y; (n) «TT. 



III. CONTINUED FRACTION 
REPRESENTATION 

In order to calculate the semiclassical transmission am- 
plitude of Eq. (||) we need to determine all families n, 

with their exiting lead (p n , q n ) and extreme points y'i . 
This evaluation naturally leads to a diophantic problem 
pof . As we show below, the coefficients (p n , q„) are part 
of the intermediate fractions (or Farey series) appearing 
in the continued fraction representation of 1/ tan 9 a ■ 

For a given angle 9 and an initial point yo, the trajec- 
tory in extended space is the straight line D(y ) defined 
by y = yo+x tan 9. Let us start with the trajectory enter- 
ing the cavity at the lowest point (yo = 0), whose exiting 
lead is defined by the segment [{po, go); {Po, go+VF)] which 
intersects D(0) (see Fig. [|). As we increase yo, we will 
remain within a family of degenerate trajectories (that 
we note by n~0) until the exiting point yo+PotB 1 n9 hits 
the uppermost point of the segment. The pair (po,qo) 
must verify the conditions: 
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a. < po tan 8 - q < W , 

b. V(p, q) such that < p tan 8 — q < po tan 8 — qo =>• 
P > Po • 

According to a, the first y-interval is [0, 5q] (or equiva- 
lently, the first y '-interval is [W— 6o, W]), while condition 
b means that (po,qo) is the first lattice point verifying a. 
The uppermost family will be associated with an inter- 
val [Oo,<?o); (Po,1a+W)], where the pair (pg.gj) verifies 
similar conditions as (po,Qo)' 

c. < W+plt&n6-ql < W , 

d. V(p, g) such that W+p% tan < IF+p tan 6-q 
W=>p>p* . 



< 



According to c, the first y-interval of the uppermost fam- 
ily is [W-8%, W] (the first y'-interval is [0,^]), while d 
implies that (Po,Qo) i s the first lattice point verifying c. 

Now that we have determined the lowest and the up- 
permost families for the [0, W\ y-interval, the following 
sequences of families can be obtained by reducing our- 
selves to the y-interval [S , W—Sq], and with the changes 
of (p , qo), (po, Qo) b y (Pi) 9i)> (PiiQi) tne conditions a-d 
define the next two families. Continuing this procedure 
until the two sequences of families begin to overlap each 
other, we obtain all the terms to be included in the sum 
of Eq. (||) . For the sequence of lower families we have 

y'iiVhQi) = W +(pi -pi_i)tan0 - (q t - , (11a) 



y' f {pi,qi)=W , (lib) 

8i=Vf- v'i = m- Qi-i - (Pi -Pi-i)taa6 , (11c) 
while for the upper-families we have 

y'M,q:)=0, (12a) 

y'f(Pl,Ql) = (P« -P«-i)tan6 - (q* u - q* u _ x ) , (12b) 



(K -p*u-i) tan61 - ill -ol-i) 



(12c) 



The very last family is simultaneously shadowe d by lower 
and u pper families, therefore has y[ given by (11a) and 
y' f by dl2b|) . 

We can now establish the relationship with the contin- 
ued fraction representation of 8 = 1/ tan 8, that is defined 
by the sequences (a m ) and (a m ) as follows |30|]: 



a a = 6 
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a = [ao] 

O-m+l — [CKm+l] 



(13) 



[a] denotes the integer part of a. The best rational ap- 
proximations to are the fractions P m /Qm, called con- 
vergents, and obtained from the recurrence relations 
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Qr 



Pm-2 - 
'■ Qm-2 



(14) 



with the initial choice of (P_i,Q_i) = (1,0) and 
(Po,Qo) = ([©])!)■ Since (P m ) and (Q m ) are sequences 
of integers, we can represent the convergents as lattice 
points (P m ,Qm) that approach the straight line D(Q) as 
m increases from above (even m) and below (odd m). 
Moreover, the convergents verify 



e. V(p, q) such that |ptan# - 

P > Pm- 



<7| < \P m taa6-Q, 



The intermediate lattice points on the segment 
[(Pm-2,Q m -2),( P m,Qm)]) define the m-th Farey se- 
quence (or intermediate fractions) by 



— Pm-2 
= Qm-2 



kP m —i 
kQ 7n —i 



< k < a n 



(15) 



The sequence of equally spaced lattice points (p m ,q m ) 
starts at the convergent (P m _2, Q m -a) (f° r fe=0) and fin- 
ishes at (P m , Qm) (for k = a m ). The translation vector is 
given by the coordinates of the convergent (P m _i, Q m —\) 
and the intermediate fractions verify the properties 

/. \fm,k,{p,q) such that < ptan# — g < 
P k 2m+\ tan 6< - g/f m+1 p > p% m+1 , 

g. Vm,k,(p,q) such that — P2m^ n ^ > Q ~ 
ptajiO > p > p\ m . 

These are the conditions b and d, respectively, de- 
termining the lattice points associated with the exiting 
leads. Therefore, the two types of families contributing 
to tb a are given by the Farey sequences of even and odd 
order convergents, with the the additional requirements 
that p n is odd and that the exiting points are closer than 
W to the straight line D(0). 

The continued fraction representation of a generic (ir- 
rational) real number results in an infinite sequence of 
convergents. However, the fact that we have defined a 
finite distance W to approach D(Q) implies that once 
the lower and upper families begin to overlap the sum in 
Eq. (||) can be cut. Therefore, for each angle 8 a the sum is 
actually finite. In addition, we show in the sequel that at 
most three sequences of intermediate fractions contribute 
to it. Let us consider the first convergent (P m , Q m ) con- 
tributing to the sum (the previous convergents being far- 
ther away than W to D(0)), and assume for definiteness 
m to be odd. Clearly, before this convergent we can only 
have one contributing Farey sequence: 



(p: 



k k 

mi y?7 



) 



with k e (1, a r , 



(16) 



The fc = case is excluded because otherwise (P m ,Qm) 
would not be the first convergent of the sum. The fol- 
lowing sequence, 



(Prn+l) Qm+l) 



with k € (1, a m +i) 



(17) 
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is relevant since it finishes at the next convergent 
(P m +i,Q m +i) = (Pm+i>9m+i)) which must contribute 
to the sum. It is closer to D(0) than (P m , Q m )) and it is 
not completely shadowed by (P m , Qm) (they are at oppo- 
site sides of -D(O)). Consequently, there are two possibili- 
ties: i) (P m +i,Qm+i) is partially shadowed by (P TO , Q m ) 
and therefore the sum stops there since the points 
(Pm+2i lm+2) with &:>0 will be completely shadowed; ii) 
(Pm+i, Qm+i) is n°t shadowed by (P m ,Qm), therefore 
(Pm+2i9m+2) is simultaneously shadowed by (P m ,Q m ) 
and (P m+ i, Q m+ i) since W+q} n+2 -p} n+2 tan9 > Q m+1 - 
P m+ itan6* > is equivalent to W + Q m — P m tan 6 > 0, 
which is precisely the condition for (P m , Q m ) being an ex- 
iting point. The points (p m+2 , lm+2) with k> 1 are com- 
pletely shadowed, and therefore the series stops there. 
In conclusion, in the case i) the contributing Farey se- 
quences are those of Eqs. (|ll)-(0), while in case ii) we 
have to add one intermediate fraction of the next family, 
namely (p m+2 ,q m+2 ). 

The above considerations simplify the semiclassical 
transmission amplitude of Eq. (||) to a finite sum, having 
at most three sequences of intermediate fractions. In 
case i) the form of the y'-intervals (Eqs. ([i"l|)-(|l2|)) is 
[W — Si,W] or [0,5*], therefore Eq. (H) can be written 
as 
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^ e n e' n exp [ikL r 



x {A n (a + b) -A„(a-b)}S„ , 



with the family-dependent function A„ defined by 



2W 
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exp 
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(18) 



(19) 



e' n = 1 if n = u (upper family) and e' n = — 1 if n= I (lower 
family). As previously stated, only odd p n should be 
considered in the expansion for the transmission ampli- 
tude. In case ii) the last family does not have a j/'-interval 
with the simple form of the previous ones. Therefore, 
Eq. (18) depending only on the weights S n can be used 
for all terms (families of trajectories), except for the last 
one, were the form (|J) together with the limits (11a) and 
(12t) should be employed. However, for the last fam- 
ily the weight S(pl n ,q^ n ) is very small and to first order 
in (a± b)S(pj n ,q^ n )/W, we can still use Eq. (||). Wc 
stress that the subindex n runs over contributing fami- 
lies, while m labels the convergents (contributing or not) 
and the superindex k orders the intermediate fractions of 
the Farey sequences. 

The forms (§h and ([l8|) of the transmission amplitude 
provide a very powerful method to numerically compute 
the conductance through a rectangular cavity. For each 
a = 1 . . . N, we only need to calculate a finite number of 
convergents and intermediate fractions of the continued 
fraction representation of 1/ tan 9 a and recursively obtain 
the weighting intervals 5 n . The advantage over quantum 



methods based on discretization (recursive Green func- 
tion, wave function matching, etc) is that it can be used 
for large wavevectors k. The method actually gets more 
exact with increasing k since it involves a semiclassical 
approximation. Obviously, Eq. ( |l8| ) docs not incorpo- 
rate diffractive effects like reflection off the lead mouths 
p|, ^0|j29| . However, as shown in Ref. p!p| , the semiclas- 
sical approach can be adapted in order to describe these 
simple diffractive effects. 
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FIG. 2. Conductance through a square cavity in units of 
e 2 /h as a function of the number of modes iV = kW/iv for an 
opening of W = 0.2 in units of the size of the square. Inset: 
local average of the conductance fluctuations (ST 2 ) versus TV. 
The straight line is the local variance of the fluctuations com- 
puted from Eq. (§]]). 

In Fig. |^ we present the conductance of a square cavity 
calculated from Eq. ( |l8| ) for an opening W = 0.2 in a large 
/c-interval (the /c-mesh shown is relatively sparse for view- 
ing purposes). We remark two salient features: the linear 
increase of the mean conductance with = kW/n and 
the fluctuations around the mean, which become larger 
as k increases. The linear increase of (g) with a slope 
given by the classical transmission coefficient constitutes 
a checking of the procedure, the increase of the fluctua- 
tions obtained within our semiclassical approach is con- 
sistent with previous quantum computations §. In the 
following section we show that further approximations on 
the semiclassical transmission amplitudes t^ a allow us to 
obtain both these features analytically. 

As it has been shown in previous studies comparing 
semiclassical and quantum results [^|, ^o|Jll| l , the former 
reproduce the gross structure of the fluctuations in the 
intermediate fc-regime (not too small k in order for semi- 
classics to apply and not too large k to prevent the quan- 
tum calculations becoming unreliable). The quantum- 
semiclassical correspondence is particularly clear at the 
level of the Fourier components of the transmission co- 
efficients. Specifically, for square cavities Wirtz et al. 
Jll| have identified the peaks of the Fourier transforms 
with the families (or bundles) of classical trajectories 
contributing in a semiclassical expansion. Unitarity, the 
mathematical translation of charge conservation, is a crit- 
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ical test for semiclassical approximations. In our case, 
the difference T + R — N grows with N, but it remains 
smaller than the fluctuations of T. 



IV. MEAN CONDUCTANCE 



The semiclassical expression (|18|) for the transmission 
amplitude is very useful since, as we have shown in the 
previous chapters, the numerical evaluation of a few con- 
tinued fractions allows the calculation of the conduc- 
tance. In the sequel we further simplify this semiclassical 
approach in order to render it analytically tractable. The 
function A„(x) defined in Eq. (1^) is peaked for x = 
(when the quantized angles of the incoming and outgo- 
ing modes coincide with the angle of the trajectory) and 
decays on the scale of W/6 n , therefore it can be approx- 
imated by the rectangular function 



n„(i) 




2W 
otherwise 



if 



x < 



W 
25~ n 



Thus, the semiclassical t ba simplifies to: 
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COS &b 

cos 9„ 



2^II n (a - b)e n e' n S n exp[ikL n ] . (20) 



The total transmission coefficient is obtained by sum- 
ming the magnitude squared of the transmission ampli- 
tudes over the propagating modes. Therefore, within our 
semiclassical approximation, it will be expressed as a sum 
over pairs of families of trajectories (with odd p n and 
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J_ y, COS^ ^ 

a,b— 1 n,n' 

x exp [ik(L n -L n >)] U n (a-b)U n r(a-b) . (21) 

From this highly oscillating function of k, we want 
to extract its secular behavior, linear in k. Averaging 
T(k)/k over all k (or T(k) over several oscillations) 
amounts to make the diagonal approximation between 
the families of trajectories (n = n'). In the well-studied 
case of isolated trajectories, the diagonal approximation 
yields the classical probability of transmission by pairing 
individual trajectories. In the present case, the concept 
of families of trajectories replaces the role of individual 
trajectories. Inserting the definition of the function II, 
we find 



(T) 



1 

W 2 



N 



EE* E 



a=l 



\b=b„ 



COS 9b 
COS 6 a 



(22) 



bmin = max{a-W/2(5„, 1} and b max = 
mm{a + W/25 n , N}. We now split the sum over fam- 
ilies according to 



(T) 



1 N 1 

— Y — 

W 2 ^—f cos ( 



N 



E 



S a <S n <8g 



E ^E cos0b_ 

=1 ~ LO<<5„<<5„ 6=1 

bmax a+W/2$„ 

cos o b + s 2 cos o b 

b=b mi „ S n >S,3 b=a-W/2S n 



(23) 



with 6 a and Sp respectively the min and max of 
{W/2a, W/2(N - a)}. In the classical limit of TV = 
kW/ir 3> 1 the last term dominates and the sum over 
b can be approximated by an integral leading to 



<d = EE 



ir 



(24) 



For each mode a we have simply obtained the total 
weight of trajectories contributing to transmission. In 
the classical limit the sum over a is converted into an 
integral over the initial angle 9 and we write 



(T) = Nt 



with 



d9 cos ( 



E8n 
W 



(25) 



(26) 



Therefore, the total transmission coefficient is propor- 
tional to the number of modes, and the constant r is 
a purely geometric factor. Breaking the contribution of 
families into that of individual trajectories we are left 
with the usual classical transmission probability 



d9 cos 9 
-tt/2 ./u 



w 



w 



f(y,o) 



(27) 



where f(y, 9) = 1 if the trajectory with initial conditions 
(y, 9) is transmitted and f(y, 9) — otherwise. 

The mean slope in the numerical results of Fig. ^ (with 
W = 0.2) is 0.67, in excellent agreement (better than 
0.5%) with the classical transmission probability r. Also, 
the local mean (T(k)) of Eq. ( p4| ) presents the same slope 
r and the mode quantization gives rise to fluctuations 
around this mean (which remain much smaller than those 
of T{k)). The classical coefficient r can be obtained, as 
in Eq. (p7[), by sampling the space of classical trajecto- 
ries with random choices of initial conditions 9 and y , 
or more efficiently, by sampling the angles 9 and incor- 
porating the weights 5 n emerging from the intermediate 
fractions of 1/ tan#, as suggested by Eq. (|26|). We then 
see that the continued fraction approach is not only useful 
for evaluating semiclassical effects, but also for classical 
properties like the transmission coefficient or the length 
distribution. Random sampling of classical trajectories 
is an appropriate procedure for chaotic structures, where 
the ergodicity of phase space results in an exponential 
distribution of lengths. On the other hand, integrable 
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cavities exhibit power-law distributions, which are more 
difficult to obtain by trajectory sampling. In this case, 
the continued fraction approach is very efficient since, for 
a given angle, only a finite number of terms are relevant, 
and the contributing families are incorporated at once 
according to their weight. 

The straight line that best approximates the transmis- 
sion coefficient T(k) is tN + k. The constant k is related 
with the elastic backscattering and also depends on the 
underlying classical dynamics. In particular, the depen- 
dence of k on the magnetic field results in the weak lo- 
calization effect |^2|,0|. However, we will not address the 
effect of a magnetic field in this work, nor the calculation 
of K. 



V. CONDUCTANCE FLUCTUATIONS 

As visible from Fig. ||, the oscillations around the mean 
transmission coefficient tN + k grow with larger N, con- 
sistently with the numerical quantum mechanical calcu- 
lations of Ref. . We will now evaluate the local fluctu- 
ations, that is (ST 2 ) = ((T - (tN + n)) 2 ) . We begin with 
the simplified expression (Elf) of T, and write 



T 



I 

W 1 



N 

E 

a.b,a'b'- 



COS Ob COSfcV >— ^ 

COS 9 a COS 9 a ' f-f. , 

-1 n.n ,n ,n 



exp 



ih(L n L n f + L n " — L r< 



x n n (a - &)ILy (a- 6)n„» (a' - b')U n ^ (a' - b') 



(28) 



If we take the average over several fluctuations, we only 
consider the terms having a null phase, i.e. 



L n — L n i + L n " Ln'" — 



(29) 



This condition is satisfied with the pairing L n = L n i and 
L n n = L n m, but the resulting term cancels against the 
square of the average transmission coefficient. A non triv- 
ial pairing is obtained when L n = L n m and L n > = L n >> 
with L n L n i, which implies a = a' . The contribution of 
this pairing to the local fluctuations is 



N N 



COS Vb COS Ob' 



W 4 ^ cos 2 Q 

a=lb,b' = l 



x Y / u n{a-b)n rl '(a-b)n n ,(a-b')Il n (a-b')Sl5 2 n , 



^EE(™{^.^}) 2 «. (30) 



In the semiclassical limit the sum over a can be con- 
verted into an integral dictating a linear behavior of 
(ST 2 )] with respect to fc 



($T 2 )i = — / docosoY 

2 J-tt/2 , 



min{(5„ , S' n } 
W 



(31) 



The proportionality coefficient is only related to the ge- 
ometry of the cavity, and can be computed using contin- 
ued fractions. 

The two pairings above described are those usually 
considered in dealing with chaotic cavities, except that 
in such cases we take individual trajectories instead fam- 
ilies. In the integrable system we are studying there is 
another non trivial way of satisfying Eq. (|29|), that is, 

L„ - L n ' = L n " - L n '" , with L n ^ L n > and L n ^ L n " ■ 
This typically happens when n,n',n", and n'" belong to 
the same (m-th) Farey sequence and they are respectively 
associated with the intermediate functions (p*j+ J ', 



Om>9m)> (Pm>2m)> and (pj£+J, g£+i), withfc ^k', j^Q 
and k, k', k+j, k'+j e (0, a m ). Also, since we are dealing 
with transmission coefficients, we need p 1 ^, pj^, p^ 3 and 
to be odd, which implies that if j is odd, P m -i must 
be even. Under such conditions, according to Eqs. ( JlO| ) 
and dll) we have 



L n - L n > = j (P m -i cos9 a + Qm-\ sin^a) 



(32) 



If neither n nor n' correspond to the first family of the 
sequence (k,k' ^ 0), Eqs. (JTTcJ), ( |T2c| ) and (|l|) dictate 
that the four exiting intervals have the same length 



\P m -i tan6> - Q m -i\ 



(33) 



This last contribution to the local fluctuations can be 
expressed as a sum over the convergents 



JY 



<*rv = £E E E^> 



a=l m k,k'=0 j 



(34) 



with the above specified restrictions for fc, k' and j. As 
in the previous case, converting the sum over a into an 
integral yields a contribution (ST 2 )u to the local fluc- 
tuations that is linear in fc, with a purely geometrical 
coefficient given by the continued fraction representation 
of 1/ tan 9. 

The linearity of (ST 2 ) = (<5T 2 )/ + (ST 2 ) H with fc that 
we have demonstrated is in good agreement with our nu- 
merical results of Fig. g. In the inset we show the numer- 
ically obtained (ST 2 ) and the straight line representing 
(ST 2 ) i from Eq. (|3l|), with a proportionality coefficient of 
0.13 calculated from the continued fractions. The agree- 
ment is very good and shows that (ST 2 )u is negligible 
in comparison with (ST 2 )j. Such behavior is understand- 
able since, for a given incoming mode a, the latter is given 
by a double sum over families of trajectories (n and n'), 
while the former is given as a [single] sum over the conver- 
gents (to) and sums within the associated Farey sequence, 
provided it verifies the above established conditions. The 
denominators Q m of the convergents increase exponen- 
tially with m, Q m > 2( m ~ 1 )/ 2 pel , or more precisely as 
Greenman has recently proved fl31| 
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lim Q 



l/m 



exp 



12 log 2 



(35) 



bar a set of measure zero. Therefore, the a m of Eqs. ( |l3| ) 
and (jlj) , giving the number of intermediate fractions of 
the m-th sequence, remain of order 1 even in the large 
to- limit. Consequently, there are much more terms con- 
tributing to (<5T 2 )/ than to {5T 2 ) H . 

In chaotic cavities the diagonal approximation is not 
enough to describe the magnitude of the conductance 
fluctuations, due to the exponential proliferation of tra- 
jectories that results in almost degenerate actions. This 
is clearly not the case in rectangular cavities, where we 
have shown the good agreement between (ST 2 )j and the 
numerical results. 



z 
< 
o 




FIG. 3. Correlation function locally normalized accord- 
ing to Eq. @ for five N intervals: 0-20 (thick solid), 
20-40 (thin solid), 40-60 (dotted), 60-80 (dashed) and 80-100 
(long-dashed). Inset: blow up of the small AiV region show- 
ing a triangular behavior at the origin for all the intervals 
considered. 

Conductance fluctuations are usually characterized 
by their amplitude (ST 2 ) and their energy-correlation 
length. In chaotic cavities, random matrix theory and 
numerical simulations [|l9|,^0) yield an universal (ST 2 ), 
independent of energy. The correlation length can be 
obtained from a semiclassical approach |l8|,^| and is pro- 
portional to the escape rate (the inverse characteristic 
time that scattering trajectories spend in the cavity). In 
our integrable cavity we have seen that (ST 2 ) is not uni- 
versal, but energy dependent. Also, there does not exist 
a characteristic time for exiting the cavity. Therefore it 
is not obvious that a correlation function depending only 
on the energy increment can be defined. That is why we 
consider the correlation function 



C(Afc) = 



(ST{k+Ak)ST{k)) 
WW) 



(36) 



where k varies on an interval much larger than Afc, but 
small enough to neglect secular variations. 

From Eqs. (|3^) and ( pf ) we obtain numerically the cor- 
relation functions shown in Fig. || for various N intervals 



from 0-20 (thick solid) to 80-100 (long-dashed). Except 
for the first interval, where the semiclassical approxima- 
tion is most questionable, the form of C(Afc) seems to be 
almost independent of the region of evaluation. More- 
over, a singularity for small Ak appears in all the N- 
intervals in the form of a cusp around the origin (in- 
set). The linear behavior of C(Afc) is to be contrasted to 
the lorentzian correlation function expected for a chaotic 
cavity. This difference between integrable and chaotic 
cavities is reminiscent of what happens with the weak 
localization effect: a semiclassical approach predicts a 
linear magneto resistance for integrable cavities and a 
lorentzian lineshape for chaotic ones |22| . The linear be- 
havior of the correlation function is consistent with the 
quantum calculations of Ref. jfjj that yielded a decay of 
the power spectrum (Fourier transform of C(Afc)) as x~ 2 
at large x. 

The oscillatory structure of C(AN) on the scale of 
AN ~ 0.2 corresponds to a typical length of twice the 
size of the square. We can understand this behavior by 
noticing that, for a given angle 9, the smallest possi- 
ble difference between the length of two trajectories is 
2. The analytical calculation of C(Ak) is considerably 
more involved than that of (ST 2 ); we will have a sum 
as in Eq. j30|), where each term should be multiplied by 
exp [iAk(L n — L n i)]. Therefore it looks rather difficult to 
extract analytically the triangular behavior around the 
origin. 



VI. CONCLUSIONS 

In this work we have developed a semiclassical method 
to calculate the transmission through a rectangular cav- 
ity, based on continued fractions. We have proved that 
the expansions, over families of trajectories, are finite 
and reproduce the basic features of quantum mechani- 
cal calculations. The conductance fluctuations, within 
our semiclassical approximation, are shown to be qual- 
itatively different from the chaotic case: they increase 
linearly with the incoming Fermi wave- vector (with a pro- 
portionality coefficient that can be easily obtained with 
our scheme) and the correlation function presents a cusp- 
like singularity at the origin. The continued fraction ap- 
proach is also useful to address classical properties like 
the transmission coefficient, the length distribution, or 
the area distribution (responsible for the the shape and 
magnetic field scale of the weak- localization effect). 

We have seen that open systems exhibit larger fluc- 
tuations when the classical dynamics is integrable. The 
situation is analogous to the density of states of closed 
systems, which is characterized by stronger fluctuations 
in integrable than in chaotic geometries. The augmented 
fluctuations in integrable closed and open geometries can 
be traced to the same origin: the bunching of trajectories 
into families in the semiclassical expansions, the Berry- 
Tabor formula and Eq. (|lj) respectively. 
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In the example discussed in this work the classical dy- 
namics in the leads has the same conserved quantities 
than in the cavity, leading to bunching of trajectories 
in families of degenerate action. This is not the case in 
circular cavities, where a semiclassical expansion based 
on individual trajectories has been developed |9). There- 
fore, the geometry of the leads plays a very important 
role Pq] and renders the quantum signatures of integra- 
bility in open systems quite involved. 

The existing experimental results 1 12 17|] show the dif- 
ference in the conductance fluctuations of integrable and 
chaotic cavities, but more work should be undertaken to 
test our results. The variation of Fermi energy (or num- 
ber of modes N) has been achieved in chaotic and in- 
tegrable ballistic microstructures jl^,^2|. However, the 
intervals AiV of variation remain much smaller than 
those allowing to detect a linear increase of the vari- 
ance. Therefore, microwave cavities |l5| appear as more 
promising for the observation of such an effect. 
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